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We examine the transport behaviour of non-interacting particles in a simple chan- 
nel billiard, at equilibrium and in the presence of an external field. The channel walls 
are constructed from straight line-segments. We observe a sensitive dependence on 
the model parameters of the transport properties, which range from sub-diffusive 
to super-diffusive regimes. In non-equilibrium, we find a transition in the transport 
behaviour between seemingly-chaotic and (quasi-) periodic behaviour. Our results 
support the view that normal transport laws do not need chaos, or quenched disor- 
der, to be realized. Furthermore, they motivate some new definitions of complexity, 
which are relevant for transport phenomena. 
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I. INTRODUCTION 

One of the fundamental aims of statistical mechanics is to shed light upon the rela- 
tionship between the macroscopic properties of a system and its underlying microscopic 
behaviour. In formalising such a relationship, conditions of "molecular chaos" are crucial, 
and are commonly assumed in order to extract physical properties from microscopic models 
of macroscopic systems (including from molecular dynamics simulations). In fact, various 
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"Chaotic Hypotheses" (CH) are either explicitly or implicitly assumed in the existing lit- 
erature; for instance Flefs.^,!^ consider the microscopic dynamics as chaotic, in the sense 
of uniform hyperbolicity, or of positive Kolmogorov- Sinai entropy, and the CH of Gallavotti 
and Cohen explicitly states: 

CH (Gallavotti-Cohen 1995). A reversible many-particle system in a stationary state 
can be regarded as a transitive Anosov system for the purpose of computing the macroscopic 
properties of the system. 

This hypothesis does not mean that particle systems actually are of the Anosov type, as the 
mathematical definition implies (^J, because the Anosov property is invoked only to obtain 
results on a limited number of quantities: those which characterize the macroscopic state 
of the system. Therefore, some of the consequences of the CH are expected to apply to 
systems which are not Anosov, the Gallavotti Cohen Fluctuation Theorem (GCFT) being 
one of these consequences jlj]. 

The Ergodic Hypothesis (EH) plays a similar role in statistical mechanics. The assump- 
tion of ergodicity in equilibrium systems allows one to make a formal association between the 
thermodynamics of a molecular system and its underlying mechanical properties, despite the 
fact that ergodicity is notoriously difficult to prove, and in many cases of physical interest 
it is known to be violated. Therefore, ergodicity is assumed for the purpose of extracting 
macrosopic data from microscopic models, largely on the basis that thermodynamic data 
thus generated appear consistent with observations of real physical systems. 

At present there is a wide body of results on the applicability of the EH, and theoretical 
explanations of why it works have been available for quite some time 0, 0] . Differently, the 
chaotic hypotheses, recently introduced to describe nonequilibrium systems, are not well 
understood yet, and the question of how much microscopic chaos, and which notion of chaos, 
is sufficient to explain the observed macroscopic behaviouir remains elusive For instance, 
recent results indicate that thermodynamic-like connections between the microscopic and 
macroscopic nature of a particle system can be forged for weaker-than-chaotic dynamical 
systems. In particular, papers such as Refs.j?! ^] show that behaviours which look like 
normal diffusion and heat conduction can be observed in polygonal billiards, although their 
dynamics are not chaotic. Indeed, their topological entropy vanishes, which means that 
nearby orbits do not separate at an exponential rate. Nevertheless, pairs of orbits eventually 



separate, with the only exception being pairs of parallel periodic orbits Consequently, 
these systems exhibit a certain sensitive dependence on the initial conditions, and their 
dynamics may appear highly disordered, indistinguishable to the eye from chaotic motions. 
Similarly, diffusive behaviour is observed in non-chaotic one dimensional lattices of maps, 

n nn 

and in the Ehrenfest wind-tree model [41], with spatial quenched disorder |11L 1121] . Further- 
more, the GCFT has been verified for an initial diffusive phase, durin g tr ansport observed 
in a nonequilibrium version of the periodic Ehrenfest wind-tree model |l3j, and in a model 
of a mechanical "pump" , with flat sides 3| • 

It is, therefore, of interest to explore such non-defocussing systems in greater detail, in 
order to understand more deeply the nature of transport in non-chaotic particle systems, and 
to determine its dependence on the geometry of the medium in which transport takes place. 
In addition to these more fundamental issues, there is also the practical consideration of the 
study of transport in porous media, as technological developments lead to the production of 
transport membranes designed ever more precisely on the atomic scale. In systems where the 
molecular size is of the same order as the pore width, solid-fluid interactions can dominate 
intermolecular interactions at laboratory temperatures and pressures, to the extent that 
the contribution from these intermolecular interactions to the overall transport rate can be 
neglected, and it is sufficient to consider only the solid- fluid interaction In this respect, 
the systems studied in this paper resemble transport in microporous membranes at low (but 
practically relevant) densities, where interactions are rare. However, such porous media lie 
at the fringes of our understanding, in terms of connecting their microscopic dynamics and 
their macroscopic properties, so that their applications can be further developed from an 
improved understanding of these theoretical issues. 

In this paper, we consider the mass transport of point particles in a simple two- 
dimensional polygonal channel, described in section |Hj consisting of a unit cell that is 
replicated infinitely in one dimension, and bounded in the other dimension. Despite be- 
ing arguably the simplest particle system that could be conceived, we observe (in section 
HVj) a surprising array of transport behaviours — while certain properties of the system 
demonstrate almost trivial behaviour, other properties display an unpredictable richness, 
which can be expressed through a sensitive dependence of the macroscopic behaviour on the 
parameters which define the geometry of the boundary. 

In our view, the dependence of the mass transport on the geometry of the system charac- 
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terizes such transport as "complex". Notions of comp 
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exity of billiards dynamics, based on 



12 1 , but in the context of the present 



symbolic dynamics, are commonly considered, cf. 
paper, it seems more appropriate to focus on the complexity of the mass transport as a 
function of the parameters of the system, rather than the complexity of the dynamics for 
one given set of parameters. Concerning this aspect, there are various studies which show 
that noninteracting particle systems, or lattices of maps, have highly irregular dependences 
of their transport coefficients on the values of the paramenters which define them (see, for 

instance, 0,0, 0,0,0)- 

These studies concern chaotic dynamical systems, and show 
that the diffusion coefficient, or the conductivity, can be quite irregular functions of the ap- 
plied field, of the shape, or of the map's slope. However, in systems of particles where these 
irregularities has been observed, their effects do not preclude important transport properties 
that are well-known in more regular systems. For instance, the irregular behaviour does not 
necessarily prevent a linear regime, close to equilibrium 23]. 

In non-chaotic systems of particles, like the ones considered here, the level of irregularity 
can grow one step further; it can produce discontinuous and, in fact, even singular transport 
coefficients as functions of the geometry of the systems. If this is the case, not only is 
the behaviour of the single trajectories unpredictable, and the dynamics is "complex" 24], 
but the overall transport itself becomes largely unpredictable. In other words, in the case 
of chaotic dynamics, the unpredictability of the fate of a single initial condition is not 
reflected in the behaviour of an ensemble of initial conditions, while here the behaviour of 
the ensemble itself — i.e. the object of interest in transport applications — is unpredictable. 
In fact, slight differences in the manufacturing of porous membranes may result in totally 
different transport properties. This leads us to the notion of "transport complexity", which 
will be discussed in section IVl 

The dependence of the mass transport on the geometry, i.e. of the statistical properties 
of the dynamics on the parameters defining the system, must be connected with the struc- 
ture of the ergodic and mixing properties of these systems. The ergodic properties are only 
known for systems corresponding to various subsets of the parameter space — however, we 
have no constructive way of explicitly specifying which parameter subsets (cf. Ref . |25| ) . 
Moreover, the ergodicity alone does not distinguish very finely between different transport 
behaviours. There are various numerical works on the ergodic, mixing and transport prop- 
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models are limited in number, the triangular billiards being the most studied models, and 
many questions concerning them remain open. Our results can be summarized as follows: 

• We find that diffusive behaviour is straightforward to obtain, despite the absence of 
chaos (in the sense of positive Lyapunov exponents), and the absence of quenched 
disorder. 

• In the nonequilibrium case we find that non-chaotic dynamics makes the use of Gaus- 
sian thermostats physically inappropriate. Indeed, the existence of a linear regime is 
not guaranteed in such systems. 

• We find a strong dependence of the transport laws on the geometry of the system. This 
shows that these non-chaotic models of noninteracting particles cannot be considered 
as proper thermodynamic systems. 

• We propose three ways of quantifying the dependence of the transport laws on the 
geometry of the systems, i.e. their transport complexity. 

II. SYSTEM DETAILS 

Let us begin by recalling the notion of a polygonal billiard. 

Definition 1. Let V be a bounded domain in the Euclidean plane M 2 or on the standard 
torus T 2 , whose boundary dV consists of a finite number of (straight) line segments. A 
polygonal billiard is a dynamical system generated by the motion of a point particle with 
constant unit speed inside V , and with elastic reflections at the boundary dV. 

As usual, elastic reflection means that the angle between the incoming velocity and the 
normal to dV, at the collision point, equals the angle between the outgoing velocity and the 
same normal. In the general theory of polygonal billiards, V is not required to be convex or 
simply connected; the boundary may contain internal walls. If the trajectory hits a corner 
of V, in general it does not have a unique continuation, and thus it normally stops there. 

In continuous time, the dynamics are represented by a flow {S ,i } tg (_ 00)00 ) in the phase space 
M = V x r 1 , where T 1 is the unit circle of the velocity angles Because the dynamics are 
Hamiltonian, the flow preserves the standard measure dxdyd'd (with x, y G M or T 2 , such 
that (x,y) EV). In discrete time, the dynamics are represented by the "billiard map" 0, on 
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the phase space $ = {(q, v) G M. : q G dA4, (v,n(q)) > 0} where is the inward normal 
vector to DM. at q, and (•, •) is the scalar product. Therefore <fr is the first return map, and 
the ^-invariant measure on $ induced by dxdydd is sin dddds, if s is the arclength on dV. 
Unfortunately, these standard measures are not necessarily selected by the dynamics, in 
the sense that different absolutely continuous measures are not necessarily evolved towards 
them. This weakens considerably the importance of the standard measures in the case of 
non-ergodic polygons. 

Definition 2. A polygon V is called rational if the angles between its sides are of the form 
■Km/n, where m,n are integers. It is called irrational otherwise. 

While polygonal billiards are easily described, their dynamics are extremely difficult to 
characterise. For instance it is known that rational polygons are not ergodic, and that they 
possess periodic orbits. But it is not known whether generic irrational polygons have any 
periodic orbit. On the other hand it is known that irrational polygons whose angles admit 



a certain superexponentially fast rational approximation are ergodic [251] . In particular the 
ergodic polygonal billiards are a dense G$ set in any compact set Q of polygons with 
a fixed number of sides, such that the rational polygons with angles with arbitrarily large 
denominators are dense in Q j^J. It is not the purpose of this paper to review exhaustively 
the properties of polygonal billiards, therefore we refer to the cited literature for further 
details. 

The special class of polygons that we consider consists of channels that are periodic in the 
x direction, but bounded by walls in the y direction, as depicted in Figure ^ The walls 
consist of straight edges, and are arranged in a saw-tooth configuration such that the top 
and bottom walls are "in phase", i.e. the peaks of the lower and upper walls have the same 
horizontal coordinate. The channel can therefore be represented as an elementary cell (EC), 
replicated along the x-axis. We denote by h the height of this cell, and set its length to 
2 Ax = 1 The heights of the isosceles triangles comprising the "teeth" along the top 
and bottom cell walls are denoted Ay t and Ay b respectively. We also introduce the mean 
interior channel height d, defined as d = h — (Ay t + A^)/2, which is equal to the mean 
height of the pore volume accessible to particles inside the channel. For convenience, we 
introduce the angles 

0i = tan- 1 , for i = b,t . (1) 
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In this paper we consider a range of values of h, Ay t , and Ay;,, which can be classified 
into two groups. First, we consider systems where Ay t = Ay^ (and consequently 9 t = 6b). 
In this case, the top and bottom saw-teeth are parallel with one another. Alternatively, we 
consider systems where Ay t and Ay& are unequal, and the saw-tooth walls are not parallel. 
The choice of Ay t = 0, Ayt ^ (or vice versa) is a special case of this last group - in 
this case, the horizontal top wall induces a vertical symmetry such that the dynamics is 
isomorphic to a system where both top and bottom saw-teeth have height Ay^, but where 
the saw-teeth are Ax out of phase, as in the equilibrium mechanical pump of Ref.[14]. We 
will therefore consider these walls as an extension of the first group of parallel saw-tooth 
walls. 

The model described so far can be called an equilibrium model, because there is no 
dissipation of energy, due to external forces or otherwise. A nonequilibrium model can be 
constructed introducing the influence of an external field e, that accelerates the particle in 
the positive x direction, and of a Gaussian thermostat, which balances the effect of the field, 
dissipating energy in such a way that the speed of the particle becomes a constant of motion. 
In this case, the particle obeys the equations of motion 



a = ep x (2) 



x = p x fix = -ap x + e 

V = Py Py = -®Py 

until it reaches the boundary dV, where it undergoes a specular reflection, as in the equi- 
librium case. The effect of the field is to curve the trajectories, making them concave in the 
direction of the field. In fact the solution of the equations of motion, for the free flight parts 
of the trajectory are given by: 

tan -A2 = tan — e~ et 3 
2 2 w 

1 , sin^(t) . A . 

x(t) = x - - In — -Al 4 
e sin v 

m - 

y{t) = yo — , (5) 

which depend on the initial angle ?V We note that the boundary of this system is not 
defocussing, and the external field has a focussing effect, so that the overall dynamics should 
not be chaotic, although it is not obvious that this is the case for all values of e. When the 
external field e is set to zero, one recovers the usual equilibrium equations of motion. 
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The data in this paper were obtained from molecular dynamics simulations of the chan- 
nel transport system described above. In both the equilibrium and nonequilibrium cases, 
the momentum and position of the particle are determined from analytic solutions of the 
equations of motion. Reported values correspond to averages of distinct simulation runs 
with different initial conditions for the particle. These initial conditions were randomly gen- 
erated with a uniform spatial distribution, and a "circle" velocity distribution, expressed by 
a probability density 

— 5{p — \) dp dd , with velocity v = p(cos sin (6) 
2n 

where the speed is \v\ = p = 1. 



III. THEORY 



Generally, our interest in studying ensemble properties stems from the foundations of sta- 
tistical mechanics. From a practical viewpoint, when studying a molecular system, we expect 
that the specific initial microscopic conditions do not alter the thermodynamic properties 
of a system. That is, we expect an equivalence between the ensemble and time averages of 
trajectory properties, such that if we were to wait long enough, the average of a property 
along almost any trajectory would be independent of the initial condition (and equal to the 
average of the property on the ensemble of initial conditions). Such a phenomenological 
requirement is incorporated into the theoretical structure of statistical mechanics through 
the mathematical notion of ergodicity, which can be summarized as follows. Consider a 
particle system constituted by N classical particles, described by the equations of motion 

x = G{x) ■ x = (q, p) G M. C M 6N , (7) 

where M. is the phase space, and the vector field G contains the forces acting on the system, 
and the particles interactions. Denote by S t x, t e M, the solution of (J7J) with initial condition 
x. The macroscopic quantity associated with an observable, i.e. with a function of phase 
$ : £1 — > R, is defined by: 

$(*) = lim - f $(5*x)* • (8) 

in which the time average represents the fact that macroscopic observations occur on time- 
scales which are very long compared to the microscopic time-scales, so that the measurement 
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amounts to a long time average. In general, however, computing that limit is not a trivial 
task at all. The problem is commonly solved by invoking the EH, which states that 

ty*) = ^J*(v)Mv) = (*)* (9) 

for a suitable measure /x (the physical measure), and for //-almost all x G Q. In principle, 
ergodic theory should identify the cases verifying the EH, and the physical measures /a, but in 
practice, this is too hard, if not impossible, to do. Nevertheless, there is now a vast literature 
on the validity of the EH, which can be understood in different ways, but which finds in 
Xhinchin's arguments on the properties of sum variables, the most convincing explanation 
^J. In practice, the time average of the functions of physical interest in systems of many 
particles is reached before a trajectory has explored the whole phase space (which would take 
too long), because such functions are almost constant, and equal to the ensemble average. 
Therefore, the finer details of the microscopic dynamics are not particularly important, 
except for the requirement of some degree of "randomness" (in order to introduce a decay 
in correlations between particles in the system). However, it is not clear what properties 
should be imposed on the dynamics in order to obtain sufficient randomness. Consequently, 
in a system devoid of dynamical chaos, such as that under investigation in this paper, it is 
of interest to investigate the behaviour of both the individual and ensemble properties of 
particles (and their trajectories). 

From a thermodynamic viewpoint, diffusion is the transport process generated by gradi- 
ents in the chemical potentials of the different chemical species in multicomponent thermody- 
namic systems, or of tagged particles moving in a host environment made of (mechanically) 



identical particles (self diffusion) 3lL l32j]. More commonly (and often more conveniently), 



this process is described in terms of Fickian diffusion, relating the mass flux to gradients in 
local density. Through the fluctuation-dissipation relation (which rests on the EH), the same 
diffusion properties are also related to the relaxation of local mass-gradient fluctuations at 
equilibrium, and therefore to the manner in which the positions of particles evolve during 
equilibrium dynamics, although the response to an external action, and the spontaneous 
equilibrium fluctuations are rather different phenomena 1441 . Fick's first law for diffusion is 

n 

expressed by [31]: 

J(x) = -DVn(x) (10) 
where J is the mass flow, D is the (Fickian) self diffusion coefficient, n is the number density, 
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and x is the position in space. This law, which can be justified in kinetic theory |32| , provides 
the phenomenological basis for the mathematics of diffusion in molecular systems, leading 
to the second-order PDE 

^ = Z^. (11) 

at dx 2 y } 

known as Fick's second law, where t is the time variable. The well-known Gaussian evolution 

n(x,t) = (4nDt)- 1/2 e- x2/4Dt (12) 

results from an initial delta-function distribution, and the linearity of (jllj) ensures that the 
diffusion of a system of molecules can be considered as the evolution of a superposition 
of Gaussians. In particular, we recover from (|12j) the linear growth in the mean-square 
displacement for macroscopic diffusion processes 

/oo 
x 2 n(x 1 t)dx = 2Dt. (13) 
-oo 

We note that, if n(x,t) is not a slowly-varying function, higher-order corrections may be 
introduced. The next approximation has the form 

dn B 2 n d 4 n _ A . 

m =D w +B 8^- (14) 

where B is called the super Burnett coefficient. In this case, the diffusion coefficient can 
still be defined as in (j!3)l — furthermore, the super Burnett coefficient can be determined 
via the relation 



{x\i)) -3<x 2 (t)) 2 = / a; 4 n(x,t)dx-3 



x 2 n(x, t)dx 



2 



2ABt. (15) 



As such, the super Burnett coefficient can be seen as a measure of the degree to which trans- 
port is diffusive, in the Fickian sense. We note, however, that it has become customary to 
call diffusive any phenomenon displaying a linear relation between the mean square displace- 
ment and the time, like (|13|). even if there are no multicomponents, or any kind of particle 
systems at hand. In this paper, we do the same. However, we note the key role played by 
the assumption of a phenomenological law, such as Fick's, in the preceding argument. In 
general, for a differing phenomenology, one cannot expect the resulting transport processes 
to remain diffusive in nature. In the absence of intermolecular interactions (or indeed, of 
other molecules), there is no phenomenological basis for expecting the mass flux to depend 
on density gradients, and it is therefore of interest to examine the resulting transport. 
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In characterising the transport law, we consider the behaviour of the displacement as a 
function of the time t, denoted s x (t), at equilibrium. In general, we expect an (asymptotic) 
relation of the form 

(sl(t)) ~ AF (16) 

where the coefficient A represents a mobility, and the exponent 7 indicates the corresponding 
transport law. The symbol (•) indicates an ensemble average, which in equilibrium systems 
is universally derived from the equal a priori probability assumption, or the EH. In the 
following we adopt the same assumption, when considering equilibrium systems, although 
it is not obvious that one should necessarily do so. Similarly, we adopt the language of 
thermodynamics, and define the following transport properties: 

Definition 3. Assume that lim^oo (s x (t)) jt 1 = A, for some A G (0, 00), then 
i. If the exponent 7 equals 1, the transport is called diffusive; 

ii. 1/7 > 1, the transport is called super- diffusive and, in particular, it is called ballistic 
2/7 = 2 (the mean square displacement is proportional to time); 



iii. The transport is called sub-diffusive if j < 1. 



Away from equilibrium, it is less straightforward to distinguish between the various transport 
laws, as the thermostat imposes an upper limit of linear growth for s x , even for super-diffusive 
transport processes. However, for the diffusive process, we expect the diffusion coefficient, 
estimated from nonequilibrium transport for a given external field of strength e, to converge 
to a well-defined value in the zero-field limit. To make this explicit, we define the quantity 

= kT{v^^ 
me 

for a system of particles of mass m at temperature T (with Boltzmann's constant k), and 
mean streaming velocity of (v x ), as the finite-field estimate, at field e, of the (nonequilibrium) 
diffusion coefficient. The existence of a linear regime guarantees the convergence of D(e), in 
the zero- field limit, to a diffusion coefficient D — that is, 

Z} = lim£>(e) = — lim^ (18) 

e^O m e^O e 

Thus, for a diffusive process, we expect to observe a linear response to the external field. 
Furthermore, in a thermodynamic system, we expect this diffusion coefficient to be equal 
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to that describing equilibrium diffusion. For a super- diffusive process, however, we do not 
expect a linear response, but rather that D(e) —>■ oo in the zero-field limit. Due to the 
thermostat, the maximum mean value of v x is equal to the initial speed of the particle v. 
For initial velocity set to unity, D(e) has upper bound of l/2e (in reduced units) for any e, 
so that -D(e) can only diverge in the zero-field limit. 

IV. RESULTS 

In this section we outline the results we have obtained from simulations of the transport 
of molecules in the saw-tooth channel system described in section |H] We examine the equi- 
librium transport properties in section IIV A[ and the nonequilibrium transport properties 
in section IIV Bl In both cases, we will be interested in the collective, ensemble behaviour 
of particles in the system, and how this ensemble behaviour relates to the behaviour of 
individual trajectories. 

A. Equilibrium 

The ergodic properties of our systems are not obvious, therefore, there seems to be no 
immediate choice for a probability distribution in phase space, to be used for the ensemble 
averages. Nevertheless, the uniform probability distribution (Lebesgue or Liouville measure, 
defined above) is invariant, and one could think that it is appropriate for transport in a 
membrane which receives particles from a reservoir, inside which the dynamics is chaotic. 
Therefore, the ensemble averages in this section are all computed assigning equal weight to 
all regions of phase space. 

1. Parallel walls, collective behaviour 

In Figure El we depict the mean-square displacement, as a function of time, for a series 
of parallel saw-tooth systems (Ay t = Ay^ = Ay), and we examine systems where the ratio 
Ay/ Ax varies from 0.25 to 3, so that the angle 8 = 8 t = 6b the saw-tooth makes with the 
horizontal varies from about 0.087T radian (rs 14°) to about 0.4-7T radians 72°). Each graph 
shows results for a single value of Ay/ Ax, for various pore heights h. For each choice of Ay 
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(recalling that Ax = 0.5), we examined pores with heights h = 1.5Ay, 2Ay, 2.05Ay, 3Ay 
and 21 Ay. The corresponding interior pore heights are d = 0.5 Ay, Ay, 1.05 Ay, 2 Ay and 
20Ay. We note that h = 2Ay corresponds to the critical pore width, above which the 
billiard horizon is infinite (i.e. there is no upper bound to the length of possible molecule 
trajectory segments without boundary collisions). We considered sets of initial conditions 
ranging from 1000 to 5000 particles. Not surprisingly, given the clearly non-diffusive nature 
of the transport observed in Figure El the corresponding super Burnett coefficients do not 
appear to have a well-defined value, but diverge over time. 

The same data have also been generated for a series of systems with one flat wall and one 
saw-tooth wall, which always have an infinite horizon. For these systems, we set Ay t = 0, 
and considered the same ratios Ay& / Ax as were examined in the parallel saw-tooth systems 
with pore heights h = 1.05 Ayt, 2.5 Ay^ and 20.5 Ay^ (which have corresponding interior 
heights d = 0.55 Ay, 2 Ay and 20 Ay). Sets of initial conditions ranged from 1000 to 5000 
particles. There appeared to be a longer initial transient period for the systems with one flat 
wall, but otherwise the results were qualitatively the same as for the parallel wall, despite 
the infinite horizon. 

Table |l] shows the values of the exponents obtained from fitting the data in Figure El and 
those for one flat wall to (|16|). for both the parallel saw-tooth systems and the systems with 
one flat wall and one parallel wall. In all cases, we observe that the transport is significantly 
super- diffusive, but not ballistic. For the Ay/ Ax = 1 system, which is a (parallel) rational 
polygonal billiard, we observe an exponent 7 ~ 1.65. All other choices of Ay/ Ax represent 
(parallel) irrational polygonal billiards. For Ay /Ax = 1/4 (where 9 < 7r/4), we observe an 
exponent of 7 ~ 1.85 in all systems. For Ay/ Ax = 2 and Ay/ Ax = 3 (where 9 > vr/4), 
we observe a similar exponent in the systems with pore height less than or in the vicinity 
of 2Ay, and a reduction in the transport exponent as the pore height increases. The same 
value 7 ~ 1.85 has been found by M. Falcioni and A.Vulpiani in the equilibrium version of 
the periodic Ehrenfest gas of Ref . [l^ . j^j]. 

We have also examined the distribution of the total x displacements s x (t), as a function 
of time. Distributions for the Ay/ Ax = 1 system (obtained from 2000 initial conditions) 
and the Ay/ Ax = 2 system (obtained from 5000 initial conditions), obtained at the end of 
the simulations, are shown in Figure El The results for the Ay/ Ax = 2 system are typical of 
the results observed for the other parallel irrational polygonal billiards. Errors are estimated 
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from the frequency counts used to generate the histograms. 

From an initial distribution that is effectively a delta function on the scale of the motion, 
the transport process produces symmetric distributions, distinct from a Gaussian distribu- 
tion. For the parallel irrational polygonal billiards, the distribution appears Gaussian out 
to one standard deviation — however, the distribution at larger displacements is clearly 
underestimated by the Gaussian. Attempted fits using functions of the form exp{— x a } fail 
to capture the shape both at the centre and in the tails, although the tails appear to be well 
modelled by a distribution of the form exp{— x L4 } (Figure^)- When 6 = 7r/4, however, the 
distribution bears little resemblance to a Gaussian (Figure Eb)- Analogous results in each 
case were obtained from the systems with one flat wall and one saw-tooth wall. 

We have also examined the behaviour of the momenta of particles in our systems. Given 
that the speed is preserved by the dynamics, the momenta can vary only in orientation, and 
we therefore examine the effect of the dynamics on the distribution of these orientations. 
Figure 0] shows the typical behaviour of the distribution of momenta orientations at the 
beginning, at the midpoint, and at the end of a typical simulation of 5000 particles. We also 
show a mean distribution obtained from averaging the momentum data over all sampled 
times, as well as over all trajectories. We find that there are no significant correlations in 
the distributions of the momenta over the course of the simulation — while the distribution 
of orientations, as a function of time, does not converge to the uniform distribution over the 
time-length we have considered (as we might expect it to for a large system of interacting 
particles), any memory effects do not appear to have a significant influence on the overall 
distribution at any instant, which remains close to uniform. Again, the same conclusion can 
be drawn from similar examination of the systems with one flat wall and one saw-tooth wall. 

Remark 1. That the angle of the walls with the horizontal plays a role in the determination 
of the transport law, while the width of the pores does not appear to do so, is reminiscent 
of the case of dispersing billiards. There, the infinite horizon adds a logarithmic correction 
to the dependence of the mean square displacement on time, which is rather hard to detect 
numerically. Despite the lack of chaos, our non- dispersing billiards could enjoy similar 
logarithmic corrections when the horizon is infinite. More surprising is the fact that in some 
cases the transport exponent (not just the matter transport) decreases as the pore width 
increases while the horizon is finite, and only increases again once the horizon is infinite. 
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2. Parallel walls, individual behaviours 

From the results we have obtained, there appear to be well-defined collective behaviours 
that are attributable to the systems we have studied — that is to say, the mean values of 
the properties obtained from simulation appear to converge, in the limit of a large number 
of independent trajectories (or particles), to well-defined values. In analogy to what we have 
presented above, we examine the evolution of the individual particle momenta, and the x 
displacements. 

The picture of the momenta is trivial in the Ay/ Ax = 1 systems, because for 9 = ir/A only 
four orientations per trajectory at most are possible, as determined by the initial condition. 
The picture of the momenta for irrational systems is less predictable. Over the course of the 
simulations, sequences of momenta, sampled at intervals o/10 4 time units, were collected. In 
Figure El we show such a sequence of momenta, sampled from a system where Ay/ Ax = 3, 
accumulated up to the 5 x 10 5 time units (Figure 2 x 10 6 time units (Figure Eh), and 10 7 
time units (Figure Eh)- The momenta are represented by symbols (circles) on the unit circle 
T 1 , while the lines indicate the sequence of the sampled momenta. It is clear from Figure El 
that, despite consisting of up to 1000 different sampled momenta, the number of distinct 
momenta visited by the particle is relatively small (of the order of 10-20). Furthermore, the 
growth of this set is gradual, and strongly correlated to the set of momenta that precede it 
in the sequence of momenta visited by the particle. The choice of 9 irrationally related to 
7r permits, in principle, the exploration of the whole unit circle of orientations. However, it 
is clear that the nature of the sequence of wall collisions limits the rate at which such an 
exploration of the unit circle can be achieved. This slow growth was observed for simulation 
times up to 10 9 time units (not shown here). 

In FigurelHlwe show the sequence of the first 1000 momenta, sampled every 10 4 time units, 
for six distinct initial conditions in a Ay/ Ax = 3 system. In each case, the available velocity 
phase space is gradually explored by the particle. We note that the rate and manner in which 
this exploration takes place (as indicated by the lines joining consecutive sample momenta) 
depends significantly, and unpredictably, on the initial condition, as is demonstrated by 
the visibly different structures generated by each. This feature is common to all parallel 
irrational polygonal billiards examined — for systems with parallel saw-tooth walls and 
systems with one flat wall and one saw-tooth wall. The flat wall in these latter systems 
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induces a vertical symmetry that is absent in the parallel systems, and appears to increase 
the range of momenta visited by particles — however, a similar degree of connectivity is 
observed in both cases. 

Finally, we turn our attention to the behaviour of the displacements s x (t) for individual 
initial conditions, as a function of time, which we expect to give further insight into the 
manner in which the momentum space is explored. In Figure we show s x (t) as a function 
of time over four different time-scales for a single particle trajectory in a system where 
Ay/ Ax = 0.25. Again, such a trajectory is typical of the results observed for the parallel 
irrational polygonal billiards examined. In Figure |H1 we show an analogous set of results for 
transport of a particle in a (rational) 6 = 7r/4 system. 

In Figure UK, the dynamics appear somewhat random on the scale of 10 6 time units, 
although on closer inspection it is clear that certain sections of the trajectory are repeated. 
Indeed, the occurrence of repeated segments is more obvious in Figure [7)d, where after 
an initial transience of 10 6 steps, the simulation appears to reach a periodic orbit, before 
changing to another orbit in the interval at about 3.5 xlO 6 , then reverting to the original 
orbit at around 8xl0 6 time units. On a longer time scale (Figure Eh), the transport appears 
to alternate between these two almost periodic orbits. However, on the scale of 10 9 time units 
(Figure |7|1), the evolution of s x takes on quite a different appearance, bearing a similarity 
with what one might observe for the random motion of a particle in a low-density gas. Such 
a resemblance suggests that, on this time-scale, there could be sufficient memory loss of 
the preceding momentum values for the trajectory to appear random. Even at this stage, 
however, the number of distinct momenta visited by the particle is still limited to the order 
to 10-100. 

In contrast, we recall that the transport behaviour in the 9 = tt/A system is restricted 
to a maximum of four distinct momenta, for each initial condition. Despite this apparently 
strong limitation, a richness in behaviour is still evident. In Figure El we observe strongly 
recurrent behaviour on all time-scales observed, although the nature of the recurrence varies 
on all observed time-scales, and is likely to continue to do so at larger and larger scales. We 
note that Figure [7|1 and Figure |Hfc represent the evolution of s x over the same time-length 
- however, while this evolution appears random in the irrational system, the motion is 
regular in the rational case. The regularity is less trivial in Figure with four small steps 
between the first two large steps, but only three small steps between the second two large 
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steps, showing that the apparent regularity doesn't make the dynamics easily predictable. 

In Figured, we show the distribution of displacements obtained by dividing the 10 9 time- 
unit trajectory for the Ay/ Ax = 0.25 system into segments of shorter time periods, in order 
to compare with the approximation to the ensemble distribution obtained from averaging 
over trajectories with independent initial conditions. It is clear from the figure that the 
two distributions are significantly different — while the ensemble distribution demonstrates 
the near-Gaussian properties observed earlier, the distribution from the single trajectory is 
significantly skewed (as could be expected from Figure |7|), and has a distribution that is 
much narrower than that of the ensemble. The correlations observed along the trajectory 
in Figure U\ lead to a distribution of displacements that is not at all characteristic of the 
ensemble, up to a time of one billion time units. 

In Figure EDd, we show the distribution of displacements obtained by dividing the 10 10 
time-unit trajectory for the Ay/ Ax = 1 system into segments, as done above. In constrast 
with the results for the Ay/ Ax = 0.25 system, the distribution obtained in this fashion 
shows excellent agreement with ensemble distribution of trajectories. 

Remark 2. The overall transport is ultimately slower in the rational systems than in the 
irrational systems. Furthermore, despite the limitations on the velocities that each single 
trajectory can take, which differ from trajectory to trajectory, the distributions of displace- 
ments exhibited by one rational trajectory equals the distribution of the ensemble, while this 
is not the case for the irrational trajectories I 

3. Unparallel walls, collective behaviours 

We have examined the transport properties of a series of unparallel saw-tooth systems, 
chosen such that the ratios Ayb/Ax lie in the vicinity of the golden ratio. In Figure ITU1 we 
show the behaviour of the mean-square displacement, and the estimated diffusion coefficient, 
as a function of time, for a series of unparallel saw-tooth systems where Ax = 0.5, Ay t /Ax = 
0.62, and Ay},/ Ax = 0.65. In analogy to the results above, we have examined these systems 
at the same range of pore heights based on the mean interior pore height d c = (Ay t + Ayb)/2 
at which the horizon becomes infinite — at d = 0.5d c , d c , 1.05<i c , 2d c and 20d c . 

Despite containing data from 10 4 independent initial conditions, the data in Figure EH 
exhibit features suggesting that they have not yet converged to a final result. These features 
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correspond to significant jumps in the mean-squared displacements (and consequently the 
finite-time estimate of the diffusion coefficient), resulting from short bursts of quasi-periodic 
behaviour (i.e. ballistic transport). We note from Figure ITUl that the effects of the bursts 
appears to grow as the horizon is opened. Similar jumps are observed in the time-evolution 
of the super Burnett coefficients, indicating that the bursts contribute to driving the system 
away from a Gaussian distribution. 

To give some sense of the size and frequency of these bursts, we show the distribution 
of displacements obtained over intervals of 10 5 time units, combining contributions from 
all 10 4 initial conditions, in Figure ^2 We observe excellent agreement with the Gaussian 
distribution, out to several standard deviations. However, in the tail of the distribution we 
find a non-negligible contribution from large-scale displacements, to which we attribute the 
behaviour of the super Burnett coefficients. These contributions correspond to the bursts 
observed in Figure EH and it is clear that a huge number of initial conditions would be 
required before the overall effect of this tail distribution could be realised by simulation. 

Table |H] shows the values of the exponents obtained from fitting the data from the ob- 
served systems to (jlfij) . and obtain behaviour that is close to diffusive. Errors in the Table 
were determined using a Marquardt-Levenberg non-linear least squares fit. The significant 
errors in the data arise from the ballistic bursts, so that a least-squares error estimate, 
more appropriate for random errors, may not be as appropriate here. However, the least- 
squares error estimates still provide useful information regarding the relative errors in the 
data obtained for the various systems. 

From Figure ITT1 it is clear that the bursts only occur for a small number of particles. We 
have found that, in each case, only 2 or 3 initial conditions are responsible for the 'significant' 
bursts — that is to say, if the contribution from these 2 or 3 particles (in 10000) is neglected, 
the resultant behaviour is diffusive within statistical error, and the fluctuations all lie within 
error about this mean diffusive behaviour. Within statistical error, one could conjecture 
that the effect of these 'significant' bursts is not sufficient to drive the behaviour away 
from diffusive behaviour, given the decay back to diffusive transport observed in Figure El 
after each burst. Clearly, however, it cannot be excluded that the effect of these bursts is 
to drive the transport at a rate somewhat faster than diffusive, either with an exponent 
slightly greater than 1, or with some slower correction, such as the hit correction for the 
Sinai billiard. We note, with the Sinai billiard in mind, that the bursts responsible for this 
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potentially super-diffusive behaviour are observed in systems with both open and closed 
horizons. 

In Figure we show the distribution of displacements after 10 6 time units, and after 
10 7 time units, from the trajectories of individual particles, again noting that the initial 
distribution is effectively a delta function since all particles begin from the same unit cell 
at s x = 0. Even after 10 6 time units, the distribution of displacements is very well fit by a 
Gaussian, consistent with our observations of a diffusive transport rate, and in contrast to 
the results for the parallel walls. 

As with the case of parallel walls, the momenta do not appear to be significantly correlated 
over the course of the simulation, and there is stronger evidence in this case of a convergence 
to a uniform distribution than in the parallel case. 

4- Unparallel walls, individual behaviours 

As with the parallel systems, we have considered the behaviour along longer individual 
trajectories, to compare with the ensemble behaviour. For unparallel walls, we find that 
there is strong agreement between the individual and ensemble behaviours, in terms of the 
distribution of both the momentum orientations and of the displacements. Sequences of 
momenta appear random along a trajectory, and the distribution of displacements obtained 
along a single trajectory demonstrates an excellent Gaussian fit, in agreement with the 
ensemble results. 

B. Nonequilibrium 

An alternative method of studying the transport behaviour is to consider the dynamics 
in the presence of an external field. For thermodynamic fluids, we expect that the transport 
coefficient determined in the linear response regime (i.e. in the zero-field limit) for the 
nonequilibrium fluid corresponds to that obtained from the equilibrium fluid properties. It 
is not clear, a priori, that such an equivalence holds for the system examined in this paper, 
so we examine the nonequilibrium transport properties, as a dynamical system of interest in 
its own right, as well as to compare its behaviour with that of the equilibrium counterpart. 
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1. Parallel 

For the non-zero-field estimates of the transport coefficient for the Ay = 2 system, with 
d = 0.5Ay, for the four different field strenghts e = 0.1, 0.01, 10~ 3 , 10~ 4 , we find that the 
qualitative behaviour of the particles is highly, and unpredictably, dependent on the field 
strength. At the highest field, e = 0.1, (Figure fTSk) all trajectories have similar qualitative 
and quantitative properties, demonstrating a "fluid-like" response to the external field - 
particles are driven in the direction of the field, and the fluctuations about a mean transport 
rate are small, since the field is strong. At e = 0.01, however, the trajectories exhibit two 
distinct transport phases — an initial phase where the particle trajectories fluctuate about a 
mean motion due to the driving field, and a second ballistic phase, where the trajectory finds 
a periodic orbit (Figure IT3b). At this field, two distinct orbits were noted — one consisting 
of 33 reflections, with period r = 9.6950889... and mean net speed Vf, ~ 0.31, the other 
consisting of 39 reflections, with period r = 12.393386... and mean net speed of Vb ~ 0.24. 
The orbits are shown in Figure El In particular, we note the existence of distinct periodic 
orbits to which the different trajectories converge, demonstrating that the dynamics at this 
field strength is not ergodic. At lower fields, a transition to periodic orbits is much rarer 
- however, bursts of almost-periodic orbits are observed, which decay after relatively short 
times to revert to the previous apparently random behaviour (Figure fTBtO. 

This transition from apparently random transport to ballistic transport has been observed 
previously in a similar nonequilibrium system with straight walls [13| . There, a transition 
time was observed between these two transport behaviours, which varied as the inverse 
square of the external field strength. Our results appear to be consistent with such a rela- 
tionship, inasmuch as the fraction of observed periodic orbits (indicating ballistic transport) 
decreases as the field strength increases. We expect that longer simulations would produce 
a larger fraction of such periodic orbits, and suggest that the reason we do not observe any 
ballistic transport in the weakest field is due to the much larger time-scale on which such a 
transition would take place. Unlike the systems in Q], where only the ensemble behaviour 
was studied, and a smooth transition from diffusive to ballistic transport was observed at fi- 
nite times, here we consider the single trajectories, and find that their individual transitions 
are sharp and occur apparently at times not limited by any upper bound. 

We show the estimates of D(e) for the applied fields e for the parallel Ay/Ax = 2 system 
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in Column 2 of Table IIHI Over such a range of fields, we would expect to observe either 
convergence to well-defined value (the diffusion coefficient), or divergent behaviour (bounded 
by the thermostat). However, we note that D(e) does not diverge in the zero- field limit, as 
we would expect for a real fluid which exhibited superdiffusive equilibrium behaviour (such 
as plug flow). 

2. Non-parallel 

In contrast to the results for the parallel systems, we find that the individual particle 
trajectories for the Ay t /Ax = 0.62, Ay^/Ax = 0.65, d = 0.5d c system display the same char- 
acteristic behaviour for the four different field strengths e = 0.1,0.01, 10~ 3 , 10 -4 examined. 
Trajectories are typified by an initial apparently random transient behaviour, followed by a 
transition to ballistic transport in a periodic orbit. The transition time also appears to vary 
inversely with the field strength. As a consequence, the fraction of trajectories observed to 
undergo a transition within 10 6 time units decreases with decreasing field, and therefore the 
lifetime of the non-ballisitic regime grows. 

For a given field strength e, each ballistic trajectory appears to have the same mean net 
speed. For each field strength considered, this mean net speed appears to be Vb ~ 0.48. We 
find that the individual trajectories converge to a single periodic orbit that is independent 
of the initial conditions, but dependent on the field strength. These orbits are shown in 
Figure H21 We note that these periodic orbits appear to be instances of a continuous family 
of periodic orbits, converging to a limiting orbit in the zero-field limit. However, this limiting 
orbit is never observed in the equilibrium trajectories, where it is no longer an attractor. 

The existence of a family of attractive periodic orbits, converging to a precise periodic 
orbit in the zero-field limit, precludes the existence of a linear regime, and is significantly 
at odds with what one would expect of an ergodic, diffusive system in nonequilibrium. 
For such systems, we would anticipate the convergence of the nonequilibrium finite-field 
estimates D(e) to the equilibrium value D in the zero- field limit. However, the limiting 
periodic orbit with a finite net speed Vj, ~ 0.48 implies that the finite- field estimate D(e) 
must diverge in the zero-field limit. 

However, as noted above, the fraction of trajectories that turn ballistic decreases with 
decreasing field, so that in the e — ► limit we expect this fraction to tend to 0. Further- 
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more, the behaviour of the trajectories before they reach a periodic orbit remains apparently 
random, and responsive to the field. Consequently, we consider the following approach to 
constructing the linear regime. At any given time, we neglect those trajectories already 
captured by a periodic orbit (ie whose transport has become ballistic), and we define a 
finite diffusion coefficient based upon the remaining trajectories. Usually, one defines the 
nonequilibrium diffusion coefficient (following the Green-Kubo result ([18))) by taking the 
limits t — > and e — > separately — here, this is no longer possible if we wish to avoid 
periodic trajectories, so the limits must be taken simultaneously. 

Therefore, given a sufficiently large ensemble of N particles, assume that the fraction 
vie; t) of particles which are still in the diffusive regime, at time t and for field strength e, 
is well-defined. Let M(N; e; t) be the subset of indices in {1, . . . , N} of these still-diffusive 
trajectories. Assume that for every 8 > there is a field e$ > such that v{e\ t) > 1 — 8 if 
e < eg and N is sufficiently large, as seems to be the case in our simulations. For simplicity, 
assume that there is a function e = e(t; 8), such that v(e(t; 8);t) = 1 — 8/2. Finally, define 
N(8) = [21/6], where [■] represents the integer part, and choose integer I so large that N(5) 
is sufficiently large for the fraction 8/2 to be sufficiently finely realized, for any 5 > 0. 

Definition 4. For every 8 > 0, distribute at random (with respect to the Lebesgue measure) 
N{8) initial conditions in the phase space. For each of them and for fixed t, consider 

A(e;t) = ^, ie{l,...,N}, 

me 

where Vi )X is the mean x-component of the velocity of the i-th trajectory (and other variables 
defined as in fli7| )J. The nonequilibrium estimate of the diffusion coefficient, if it exists, is 
defined by 

D ne = lim lim —-— , , — - — - DAe(t\ 8): t) 

V ' V V ' h > ieM(N(S);e;t) 

_fcT 1 1 ^ li9j 

"SS me (t;8) N(6)u(e(t;6);t) t f-f Sxi[th 

K ' ' v ' K K ' h ' ieM(N(S);e;t) 

Thus, for a given choice of 8, and t, we evaluate the estimate of the diffusion coefficient 
(([17))) considering only those trajectories which have not been captured by a periodic orbit. 
Then, the limit t — > oo includes also the limit e — f 0, and the limit 8 — > includes the limit 
iV — > oo. We can now define the following: 
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Definition 5. If the equilibrium system has a finite diffusion coefficient D, and D ne = D, 
the system is said to have a linear regime. 

In Table ITTT1 we report the 'corrected' estimates D(e) for the various fields, where we 
neglect contributions from ballistic trajectories. In the third column we report estimates 
for the parallel Ay /Ax = 2 system, and in the fifth column we report estimates for the 
unparallel Ay t /Ax = 0.62, Ay^/Ax = 0.65, d = 0.5d c system. We note that the 'corrected' 
data for the parallel Ay /Ax = 2 system are consistent with a divergent trend (taking 
into account the large statistical error in the weakest field). For the unparallel system, the 
equilibrium and nonequilibrium transport properties become more consistent, although they 
are not conclusively so, because of the large error bars produced by the statistical analysis at 
low field. Such noise is typical of NEMD simulations at low field, where the signal-to-noise 
ratio becomes low. Typical NEMD simulations, however, would still converge to yield the 
same transport coefficient in the long-time limit — in our non-chaotic systems, the transition 
from random to ballistic behaviour effectively places an upper bound on the time during 
which random behaviour can be observed. The estimated errors in the diffusion coefficient 
obtained will therefore depend on the rate at which the mean transition time increases with 
decreasing field. 

Remark 3. This subsection leads us to conclude that the use of thermostats needs some 
form of chaos, whether permanent or transient, in order for a linear regime to be observed. 
This is a cause for concern in the case of thermostatted non-interacting particle systems. 

V. TRANSPORT COMPLEXITY 

The properties of the dynamics in these systems are strongly dependent on the values 
of 9 t and 6^, and their relative properties — whether they divide one another, and whether 
they divide tt. This sensitivity to the boundary conditions is strongly "unthermodynamic" in 
nature, and owes its origin to the absence of inter-particle interactions. In real molecular dy- 
namical systems, chaotic behaviour acts as a "double-edged sword" — on the one hand, the 
sensitivity to initial conditions ensures that trajectories diverge exponentially; on the other 
hand, under such dynamics almost all trajectories explore the phase space in such a way that 
they display the same macroscopic behaviour. In the current systems, the dynamics appears 
to display a sensitivity to initial conditions despite the absence of dynamical chaos. How- 
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ever, the dynamics is not always sufficiently mixing to cause all trajectories to demonstrate 
the same average behaviour. By a judicious choice of saw-tooth angles, one can produce 
dynamics which appear consistent with macroscopic notions of diffusion. Clearly, however, 
in the absense of intermolecular interactions, one can only expect to describe those aspects 
of thermodynamics which are based on purely convective processes. Such models lack the 
essential element of momentum or energy exchange that would be required to describe pro- 
cesses that involve conduction, such as heat transfer or viscous effects, and consequently one 
cannot draw conclusions on such processes based on the systems of noninteracting particles, 
such as those examined here. 

Similar considerations appear to apply to chaotic noninteracting particle systems, as 
they also have been found to have thermodynamic properties irregularly dependent on the 



geometry of the medium 



18, 



19, 



20, 



221 ] . However, in these systems, the presence of chaos 



induces a higher stability, which is manifested in the fact that the transport coefficient, but 
not the transport law, is affected by the geometry The transport coefficient, in particular, 
turns out to be a close-to-continuous function |4j| , which may, or may not, be a fractal 
curve, depending on the case. As conceptually interesting as this question is, it does not 
bear serious consequences for the predicted transport, when the parameters determining the 
geometry are only known with a finite precision. Indeed, this implies a small error in the 
transport coefficient, which results in a small uncertainty on the transported quantities. In 
other words, transport is "predictable" in these systems, hence not "complex", although 
the single trajectories may display arbitrarily high complexity (in symbolic dynamics, and 
information theoretic terms, for instance). 

Systems with an unpredictability which affects not only the single trajectories, for fixed 
geometry, or the transport coefficients, for uncertain geometry, but affects the transport 
laws themselves, would seem to belong to a class with a higher level of "complexity". The 
polygonal channel studied in this paper is such a system. Other map systems with such 
dependency have also been observed - see Ref.Qfwhece the geometry unpredictably gave 
periodic, diffusive or ballistic behaviour) and Ref.[35]. Our numerical results indicate that 
the transport law can change from diffusive to sub-diffusive regimes, or alternatively quite 
high superdiffusive regimes, after small changes of the parameters defining the geometry. To 
quantify this kind of complexity, we propose the following definitions. 
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Definition 6. Consider a transport model, whose geometry is determined by the parameter 
y, which ranges in the interval [0, h], and such that its transport law is given by 

l im iflM = A , < A < oo (20) 

t^oo fy 

with 7 a function of y varying in [0,2], when y spans [0,h]. Let Aj(y m ,yM) £ [0, oo] 
be the difference between the largest and the smallest value of 7, for y in the subinterval 
{VmiVM) C [0, h], where Aj(y m ,y M ) =00 if in (y m ,yM) there are points for which H?J) is 
not satisfied by any 7 > 0. 

i. The transport complexity of first kind of the transport model in (y m , yu) is the number 

Ci{y m ,yM) = —, -ep,oo) (21 

2{y M - Vm) 

if it exists. 

ii. The transport complexity of second kind of the transport model for y = y is the 
exponent C 2 = ^(y), if it exists, for which the limit 

C 1 (y-e 1 y + e) 

is finite. 

iii. The transport complexity of third kind of the transport model for y = y is the limit 

C 3 {y) = limA 7 (y + (23) 

These definitions are motivated by the following considerations. If C\ does not vanish, the 
system is surely highly unpredictable from the point of view of transport, in the interval 
(Vm, Vm), because its transport law is only known with a given uncertainty. This could 
be the case, for instance, of a batch of microporous membranes with flat pore walls, whose 
orientation is obtained with a certain tolerance (transport complexity of first kind). However, 
C\ may diverge around some point of [0, h], as our data seem to indicate, giving rise to an 
even higher level of transport complexity. Assuming that this divergence has the form of a 
power law, we take the power as a measure of this second kind of complexity. But even this 
level of complexity seems to be insufficient for our models, which indicate that A^^y^yM) 
could be discontinuous. For this reason, we introduce the third kind of transport complexity. 
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In order to further investigate these notions, and how they relate to our systems, we 
have considered the transport behaviour in a narrow range about two principal systems, 
taken from the examples of super-diffusive transport observed in Section IIV Al The first 
principal system is the rational parallel Ay t /Ax = 1, where we consider the transport 
behaviour in the limit that Ayb/Ax — > 1, Ayt/Ax > 1. The second principal is the irrational 
parallel Ay t /Ax = 2 system, where we consider the transport behaviour in the limit that 
Ay^/Ax — > 2,Ayb/Ax > 2. In Table [TVl we show the exponent 7 of the transport laws (as 
per (J2I11)) corresponding to the various choices for Ay t and Ay b . 

We note that, in both cases, there appears to be a strong discontinuity at Ay t Ayb - 
when the walls are not parallel, the behaviour is no longer super diffusive. In the rational 
case, the transport coefficient appears to depend unpredictably on the Ay t , but is always 
distinctly sub- diffusive. By contrast, in the irrational case, the transport appears to be 
essentially diffusive when the walls are not parallel. This behaviour is consistent with that 
observed in Section \l V A 31 for the unparallel walls, which are also irrational. 

Due to the discontinuity, the transport complexities C\ and C2 diverge in these cases, 
reflecting the high degree of unpredictability that has been observed in the previous sections, 
which can only be quantified by C3. 

VI. CONCLUSION 

In this paper, we have examined the transport properties of a dynamical system of re- 
markable simplicity — a two-dimensional channel with straight-edged walls, populated by 
non-interacting, point-like particles. Importantly, our choice of system ensures the absence 
of dynamical chaos. However, this system displays a rich variety of transport behaviours, 
belying the simplicity of the underlying dynamics. 

For the parallel equilibrium systems where the pore angles 9t and 8b are equal (8 t = 
Of, = 8), we observe super-diffusive (but not ballistic) behaviour. For the parallel irrational 
polygonal billiards, we find that the root-mean-square displacement grows approximately 
as (s^(t)) ~ t 1 ' 85 in the equilibrium dynamics. However, for the rational counterpart, 
( s x(t)) ~ ^ L65 & t equilibrium, apparently constrained by the limited range of accessible 
momenta. These results hold also for the case where 8 t = 0, 8b = 8, which are isomorphic 
to systems where the top and bottom pore angles are equal, but the saw-teeth are Ax out 
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of phase. Interestingly, these results are not strongly sensitive to the opening of the horizon 
(particularly in the case of a shallow saw-tooth angle, 9 < 7r/4). It is perhaps surprising 
that no ballistic motion was observed in these systems, given their simple structure — while 
the transport is strongly super- diffusive, it is also clearly sub-ballistic in nature. 

The behaviour for individual trajectories in these parallel-walled systems are strongly 
dependent on the initial conditions, in part because of the nature in which the dynamics 
permits (or rather, restricts) the exploration of the momentum space. This can be under- 
stood through the strong correlation of any particular value in the sequence of a particle's 
momenta with those momenta preceding (and indeed, following) it, developing trajectory 
sequences through the creation of seemingly quasi-periodic "building blocks" on various 
length scales. Despite the limitations imposed by this strong correlation, the behaviour 
of the trajectory over a particular order of magnitude in time is seemingly unpredictable 
from its behaviour over shorter time-scales, and it appears that the correlation length is 
of the order of 10 9 time units, or more, when 6 is irrational with respect to tt, and orders 
of magnitude longer when 6 is rationally related to it. To give some physical sense to this 
observation, we consider an analogy to the transport of a light symmetric gas molecule (such 
as methane) along a nanopore of width 1 nm. At room temperature, the unit velocity in our 
system corresponds to a root-mean-square velocity for methane of approximately 400 m/s, 
implying that our time units correspond roughly to units of picoseconds. Consequently, the 
correlation length in the irrational systems is at least of the order of microseconds, during 
which time the particle will have travelled of the order of millimetres. For the rational 
systems, the lower bound for the correlation time is closer to the order of seconds, although 
the mean displacement over this time is of similar order (because of the slower transport, 
reflected by the smaller value of 7). We note that the mean time between collisions would 
be much less than this, for all but the most rarified of conditions. 

For the nonequilibrium parallel systems, we observe similarly strong dependence on the 
initial conditions of the trajectories — indeed, at some field strengths we observe clearly 
non-ergodic behaviour, with the existence of distinct attracting periodic trajectories (with 
different net speeds) for a single dynamical system. It appears that the existence of periodic 
orbits is determined by the field strength, although bursts of quasi-periodic behaviour are 
observed over short time intervals. There is no clear relation between the onset time for 
ballistic (periodic) transport and the field strength. More fundamentally, there does not 
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appear to be any clear connection between the non-zero-field estimates of the diffusion 
coefficient and the zero-field (equilibrium) value (which is infinite in the case of parallel 
walls). Certainly, there is no evident divergence of D(e) in the zero-field limit, as one would 
expect for a super-diffusive thermodynamic system (such as plug flow). Although, these 
systems could have a diffusion coefficient defined as in section IV, the existence of a linear 
regime has not been clearly estabilished, because the equality between D and D ne is not 
guaranteed by our calculations. 

For the unparallel systems, where 9 t ^ 9b, we observed diffusive transport in the equi- 
librium system — (s^(t)) ~ t, with a Gaussian distribution of displacements s x consistent 
with the behaviour predicted by (|12jl. The trajectories were more sensitive to the opening of 
the horizon than for the parallel systems, with occasional short-lived bursts of much faster 
transport. In contrast to the unexpected richness in behaviour of the parallel systems, in the 
unparallel systems, the properties of individual trajectories consistently matched ensemble 
properties. This was also true of the nonequilibrium results, where there was a much clearer 
field-dependence of the transport behaviour, and in particular the transition to a periodic 
orbit. In the unparallel systems, the periodic orbit was unique for a given field (and ap- 
peared to be continuously dependent on the field), and the mean onset time of the periodic 
orbit grew with decreasing field strength. This is consistent with the results of Ref. 
which, however, were obtained for a system with parallel boundary walls. However, the wall 
angles were chosen such that tan# was close to the golden ratio (Vb + l)/2 = 1.618 . . ., so as 
to "maximise" the irrationality of the relationship of 9 to n. It is possible that such a choice 
leads to much smaller correlation times in the momentum sequence for particle trajectories, 
and consequently faster convergence to the longer time-scale behaviour more reminiscent of 
intermolecular collisions, observed in Figure [7| 

Despite the apparently ergodic behaviour, however, our results are not conclusive regard- 
ing the connection between the equilibrium and nonequilibrium transport behaviours: the 
existence of a linear response regime, where the non- zero- field estimates D(e) converge to 
the equilibrium diffusion coefficient in the zero-field limit is not clear. 

Among the systems we have observed, we find that: 

• systems with irrational walls are super-diffusive if parallel, but diffusive if rationally 
related; 
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• systems with one rational wall and one irrational wall are sub-diffusive; 

• systems with two rational walls are super- diffusive. 

Furthermore, we identify a means of quantitatively estimating the uncertainty related to 
the system transport, through the parameters used to define the geometry. Such notions 
of transport complexity are useful in distinguishing between various transport systems, the 
range of transport properties that they exhibit, and in particular the predictability of these 
transport properties. The various transport properties observed in the systems examined in 
this paper, and their strongly unpredictable nature, are reflected in the divergent values for 
the proposed transport complexities C\ and C2, and in the finite value of C3. We note that, 
for the sake of simplicity, we have focussed on the dependence of the transport law on just 
one paramenter, but our investigation reveals that transport in our models may depend in a 
counterintuitive and irregular fashion on other parameters as well, such as h. Striking is the 
case in which a growth of h implies a reduction of the transport exponent. More complicated 
notions of transport complexity may then be envisaged, but the conceptual picture outlined 
here would not change. 

This analysis leads us to the following observation. As is well known, the thermodynam- 
ical properties of a macroscopic system are not a function of the boundary of the medium; 
the properties of a fluid, such as its viscosity, do not depend on the geometry of its container, 
and, in general, the nature of transport is essentially independent of the geometry. Fermi 
expressed this concept, at the beginning of his book on thermodynamics [2a], in which he 
states: "The geometry of our system is obviously characterized not only by its volume, but 
also its shape. However, most thermodynamical properties are largely independent of the 
shape, and, therefore, the volume is the only geometrical datum that is ordinarily given. It is 
only in the cases for which the ratio of surface to volume is very large (for example, a finely 
grained substance) that the surface must be considered. This is clearly understood in the 
terms of kinetic theory, which leads to explicit expressions for the mutual or self diffusion 
coefficients in terms of only the mean-free paths A for collisions among particles, without any 
reference to the shape of the containers. Only in the case that A is of the same order of the 
characteristic lengths of the container, does this play a role; but in that case, the standard 
laws of thermodynamics cease to hold, and are replaced by those of highly rarefied gases, or 



Knudsen gases |3ll 



Because our systems do not have this extreme ratio of surface 
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to volume, our results indicate that our systems cannot be considered as thermodynamic 
systems, much as they remain highly interesting from both theoretical and technological 
standpoints. Indeed, as has been noted elsewhere [3, Q]' much care must be taken in 
the thermodynamic interpretation of the collective behaviour of systems of non-interacting 
particles. The arguments of Refs. mainly referred to nonequilibrium systems: here 

we provide some examples which show that they apply to equilibrium systems as well. All 
this can be summarized stating that a certain degree of chaos, or of randomness is required 
in the microscopic dynamics of particle systems, for them to look like thermodynamic sys- 
tems, but there are two ways in which this can be achieved. If the randomness is intrinsic 
to the fluid, the geometry of the system is a secondary issue, and the fluid behaves as a 
proper thermodynamic system. If the randomness is produced by the geometry of the outer 
environment, the geometry plays an obviously important role, and systems such as ours do 
not behave like proper thermodynamic systems. 
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FIG. 1: Symmetric saw-tooth pore model used throughout this paper. The total pore height is 
denoted h, the width of the repeated unit cell is 2 Ax = 1. The "tooth" heights on the top and 
bottom are, denoted Ay t and Ayt, respectively. 
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FIG. 2: Evolution of mean-squared displacement for parallel saw-tooth systems, for (a) Ay/ Ax = 
0.25, (b) 1, (c) 2, and (d) 3. The data are obtained from 1000-5000 initial conditions, and the 
average number of collisions is of the same order as the total time. 
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FIG. 3: Distribution of displacements after 10 6 time units for parallel saw-tooth systems, where 
(a) Ay/ Ax = 2, from 2000 initial conditions, and (b) Ay/ Ax = 1, from 5000 initial conditions. 
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FIG. 4: Distribution of momenta orientations (a) at t = 0, (b) at t = 5 x 10 5 , (c) at t = 10 6 , and 
(d) averaged over all times, for the Ay/ Ax = 1 system with d = 0.5 Ay. The dotted line indicates 
the uniform distribution. Errors are estimated from the frequency counts used to generate the 
histograms (and are hence smaller in (d), where the data is comprised of more samples). 
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FIG. 5: Momentum progression for a single trajectory up to (a) t = 5 x 10 5 (50 samples), (b) 
t = 2 x 10 6 (200 samples), and (c) t = 10 7 (1000 samples), for the Ay/ Ax = 3 system with 
d = 2Ay. 
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FIG. 6: Sequences of sampled momenta for six different initial conditions, up to 10 time units 
(1000 samples), for the Ay/Ax = 3 system with d = 2Ay. 
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FIG. 7: Example of the time evolution of a particle displacement for the Ay/ Ax = 0.25 system 
(d = 0.5 Ay), for a sequence of four different time-scales. Boxes indicate the position of the previous 
graph in the sequence. 
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FIG. 8: Example of the time evolution of a particle displacement for the Ay/ Ax = 1 system 
(d = 2 Ay), for a sequence of four different time-scales. Boxes indicate the position of the previous 
graph in the sequence. 
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FIG. 9: Histogram of displacements collected along a single trajectory, and from the ensemble of 
trajectories, for the (a) Ay/ Ax = 0.25 (d = 0.5Ay), and (b) Ay/ Ax = 1 (d = 2.0Ay) systems. 
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FIG. 10: Evolution of (a) the mean-squared displacement, and (b) the diffusion coefficient, for 
transport in systems with unparallel saw-tooth walls — Ay t /Ax = 0.62, Ay^/Ax = 0.65, d = 
0.5d c , d c , 1.05d c , 2d c and 20d c . 
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FIG. 11: Displacements over intervals of 10 5 time units, combining contributions from all 10 4 
initial conditions, for transport in the system Ay t /Ax = 0.62, Ay^/Ax = 0.65, for d = 
0.5d c , d c , 1.05d c , 2d c and 20d c : (a) within three standard deviations; and (b) over the full range 
of displacements (shown using a log- linear plot). 
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FIG. 12: Distribution of final displacements for unparallel saw-tooth systems where Ay t /Ax 
0.62, Ay b /Ax = 0.65, d = 0.54- 
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FIG. 13: Displacement versus time for ten trajectories for the parallel Ay/ Ax = 2 system, for (a) 
e = 0.1, (b) e = 0.01, (c) e = 10~ 3 , and (d) e = 10~ 4 . 
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FIG. 14: Two distinct periodic trajectories of the finite field transport for the parallel Ay /Ax = 2 
system for e = 0.01. The left hand side orbit has period r = 12.393386... time units with 39 
reflections per orbit. The right hand side orbit has period r = 9.6950889... time units with 33 
reflections per orbit. 
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FIG. 15: (Attractive) periodic trajectories of the finite field transport for the Ay t /Ax = 
0.62,Ayt/Ax = 0.65, d = d c system, for 0.1 < e < 10~ 4 . The trajectories for the two weakest 
fields are indistinguishable on the scale of the figure. Cloud shown for aesthetic reasons. 
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TABLE I: Equilibrium transport exponents: For saw-tooth boundary base triangles with height- 
width ratio Ay/ Ax. For each (mean) pore height tested, the observed exponent out to 10 6 time 
units (of the order of 10 6 -10 7 collisions) is given. The number of initial conditions used to compute 
averages ranges from 1000 to 5000, and the errors are estimated by ±0.03 in all cases. 
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TABLE II: Equilibrium transport exponents: For unparallel saw-tooth boundary base triangles 
with height-width ratio Ay/ Ax. For each (mean) pore height tested, the observed exponent out to 
10 6 time units (of the order of 10 6 -10 7 collisions) is given. The number of initial conditions used to 
compute averages varies from 2000 up to 10 4 . Numbers in brackets correspond to error estimates 
from Marquardt-Levenberg least-squares fits. 
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equilibrium 



co 



CO 



0.386 ±0.05 



0.386 ±0.05 



TABLE III: Best estimates of the finite field diffusion coefficient D(e) for 0.1 < e < 10~ 4 for the 
parallel Ay/ Ax = 2 system, and the unparallel Ayt/Ax = 0.62, Ay^j Ax = 0.65, d = d c system, 
showing the raw data and data corrected for periodic orbits. Also included, for comparison, are the 
equilibrium results for the parallel case (exhibiting super-diffusive behaviour) and the unparallel 
case. Data was collected over 10 6 time units. 
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TABLE IV: Transport law exponent 7, as identified by (|20[). observed for various choices of Ayt/Ax 
and Ayiy/Ax. There appears to be no trend toward super-diffusive transport, arbitrarily close to 
(but away from) the parallel case. Numbers in brackets indicate errors in the last figure. 



